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ABSTRACT 

Results from a large set of hydrodynamical smoothcd-particle-hydrodynamic 
(SPH) simulations of galaxy clusters in a flat ACDM cosmology arc used to in- 
vestigate the metal enrichment and heating of the intracluster medium (ICM). 
The physical modeling of the gas includes radiative cooling, star formation, 
energy feedback and metal enrichment that follow from the explosions of su- 
pernovae of type II and la. The metallicity dependence of the cooling function 
is also taken into account. The gas is metal enriched from star particles ac- 
cording to the SPH prescriptions. The simulations have been performed to 
study the dependence of final metal abundances and heating of the ICM on 
the numerical resolution and the model parameters. 

For a fiducial set of model prescriptions the results indicate radial iron 
profiles in broad agreement with observations; global iron abundances are also 
consistent with data. It is found that the iron distribution in the intracluster 
medium is critically dependent on the shape of the metal deposition profile. 
At large radii the radial iron abundance profiles in the simulations are steeper 
than those in the data, suggesting a dynamical evolution of simulated clusters 
different from those observed. For low temperature clusters simulations yield 
iron abundances below the allowed observational range, unless it is introduced 
a minimum diffusion length of metals in the ICM. The simulated emission- 
weighted radial temperature profiles are in good agreement with data for 
cooling flow clusters, but at very small distances from the cluster centres 
(~ 2% of the virial radii) the temperatures are a factor two higher than the 
measured spectral values. 

The luminosity-temperature relation is in excellent agreement with the 
data, cool clusters {Tx ~ IkeV) have a core excess entropy of ~ 200keVcm? 
and their X-ray properties are unaffected by the amount of feedback energy 
that has heated the ICM. The findings support the model proposed recently 
by Bryan, where the cluster X-ray properties are determined by radiative 
cooling. The fraction of hot gas fg at the virial radius increases with Tx and 
the distribution obtained from the simulated cluster sample is consistent with 
the observational ranges. 

Key words: hydrodynamics - methods : simulations - cluster: evolution - 
intergalactic medium - metallicity. 

1 INTRODUCTION 

Galaxy clusters are the largest virialized structures known in the universe and are consid- 
ered useful probes to constrain current cosmological theories of structure formation. X-ray 
observations of galaxy clusters show that most of the baryonic cluster mass is in the form 
of a hot ionized intracluster medium at temperatures ~ 10^ — 10* °K, with the bulk of the 
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emission in the X-ray band from bremsstrahlung processes (Sarazin 1986). The dependence 
of the X-ray emission on the square of the gas density allows us to construct cluster samples 
without the biases which may arise in the optical band. The final physical state of the ICM 
is primarily determined by the gravitational processes, which have driven the dynamical evo- 
lution of the gas and dark matter mass components during the cluster collapse. Under the 
assumption of hydrostatic equilibrium the ICM gas distribution can be modeled to connect 
the gas temperature Tx to the cluster virial mass or to the bolometric X-ray luminosity Lx- 
The X-ray temperature and luminosity functions are then predicted for a given cosmolog- 
ical model according to the standard theoretical Press-Schecter (1974) mass function. The 
observed evolutionary history of these functions (Edge et al. 1990; Henry & Arnoud 1991 ; 
Henry 1997) can then be used to put severe constraints on the allowed cosmological models 
(Henry & Arnoud 1991 ; White, Efstathiou & Frenk 1993; Eke, Cole & Frenk 1996; Bahcall 
& Fan 1998; Kitayama, Sasaki & Suto 1998). 

If the ICM has been shock heated solely by gravitational processes, the cluster scaling 
relations are predicted to obey a self-similar behavior. For instance, the Lx — Tx relation 
should scale as Lx oc T|- (Kaiser 1986). For clusters with Tx ~ IkeV there is a wide ob- 
servational evidence (David et al. 1993; Allen & Fabian 1998; Markevitch 1998) that the 
observed bolometric X-ray luminosity scales with temperature with a slope steeper than ex- 
pected {Lx oc Tx)- This implies that low temperature clusters have central densities lower 
than those predicted by the self-similar scaling relations (Bower 1997; Ponman, Cannon 
& Navarro 1999; Lloyd-Davies, Ponman & Cannon 2000). This break of self-similarity is 
usually taken as a strong evidence that non-gravitational heating of the ICM has played an 
important role in the ICM evolution (Evrard & Henry 1991; Kaiser 1991; White 1991), at 
least for cool clusters. 

The most considered source of energy which has been considered as a heating mechanism 
for the ICM is supernovae (SNe) driven-winds, which inject energy into the ICM through 
SNe explosions of type la and II (White 1991; Loewenstein & Mushotzky 1996). The energy 
input required to heat the gas at the entropy level necessary to reproduce the observed 
departure from self-similarity is estimated to lie in the range ~ 0.5 — "ikeV per particle 
(Balogh, Babul & Patton 1999; Wu, Fabian & Nulsen 2000; Tozzi & Norman 2001). Several 
authors have argued that it is unlikely that SNe can provide the required energy to heat 
the gas and have suggested active galactic nuclei (Valageas & Silk 1999; Wu, Fabian & 
Nulsen 2000), as the extra energy source for ICM heating. Bryan (2000) has proposed the 
alternative view that radiative cooling and the subsequent galaxy formation (Pearce et al. 
2000), can explain the observed Lx — Tx relation because of the removal of low-entropy gas 
at the cluster cores. 

The estimates of the energy injected from SNe into the ICM are based on the measured 
abundances of metals in the ICM. The elements that have enriched the ICM have been 
synthesized in the SNe explosions of the stellar population of the cluster. Two enrichment 
mechanisms have been proposed: supernova-driven galactic winds from early-type galaxies 
(Mushotzky & Lowenstein 1997), or ram pressure stripping of the enriched gas from galaxies 
(Gunn & Gott 1972 ). Another possibility is that of a significant contribution to the metal 
enrichment and heating of the ICM, associated with an early-epoch generation of massive 
Population HI stars (Loewenstein 2001). 

A large set of observations confirms that the ICM of galaxy clusters is rich in metals 
(Arnaud et al. 1992; Mushotzky & Lowenstein 1997; Fukazawa et al. 1998; Dupke & White 
2000a; Finoguenov, David & Ponman 2000; Fukazawa et al. 2000; Matsumoto et al. 2000). 
Analyses of X-ray spectra show that the abundance of heavy elements in the ICM is nearly 
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~ 1/3 solar. These measurements provide a strong support for the SN scenario as a heating 
source for the ICM. 

Abundance gradients have also been measured (Ezawa et al. 1997; Dupke & White 
2000b; White 2000; Finoguenov, David & Ponman 2000; De Grandi & Molendi 2001; Irwin 
& Bregman 2001). For instance Ezawa et al. (1997) found a decline in the iron abundance of 
AWM7 from ~ 0.5 solar in the centre to ~ 0.2 solar at a distance of ^ SOO-R'pc. Measurements 
of the relative abundance of the heavy elements can be used to constrain the enrichment 
mechanism of the ICM and the energy input from SNe. Analysis of the spatial distribution 
of metallicity gradients is also important to discriminate among the proposed enrichment 
scenarios (Dupke & White 2000b). From the observed metallicities one can estimate the 
SN energy released given a shape for the initial-mass-function (IMF) and a nucleosynthesis 
yield model for the SN explosions. Finoguenov, Arnaud & David (2001) have obtained the 
Si and Fe abundances using the X-ray data of a selected sample of X-ray clusters. From 
the measured abundances, they find a significant contribution to the energy per particle 
associated with SNe explosions at temperatures ~ 3keV. 

These kinds of estimates of the SN energy input suffer from the theoretical uncertainties 
in the yield models and in the form of the IMF (Gibson, Loewenstein & Mushotzky 1997). 
Moreover, these estimates can vary widely with the assumed spatial distribution of the ICM 
metals and the transfer efficiency of the kinetic energy released in a SN explosion to the 
ambient gas (Kravtsov & Yepes 2000). Hydrodynamical simulations of cluster evolution have 
the advantage over analytical methods that they take into account the dynamical evolution 
of the gas. From the results of numerical simulations of galaxy formation, Kravtsov & Yepes 
(2000) conclude that it is unlikely that SNe can provide the required energy input, even 
assuming the existence of radial metallicity gradients in the ICM. 

In order to implement self-consistently the metal enrichment of the ICM in hydrodynam- 
ical simulations of galaxy clusters, one first must consider the effects of non-gravitational 
processes on the cluster gas distribution. The effects of radiative cooling in simulated clus- 
ters have been investigated by a number of authors (Anninos & Norman 1996; Lewis et al. 
2000; Pearce et al. 2000; Yoshikawa, Jing & Suto 2000; Valdarnini 2002). One of the main 
conclusions of these simulations is that the modeling of radiative processes for the gas can- 
not be decoupled from a prescription for turning cold, dense gas into stars. This is done in 
order to avoid unphysically high densities (cooling catastrophe). The luminosities Lx of the 
simulated clusters are found to be physically plausible, provided that a suitable prescription 
for the treatment of the cold gas has been taken into account. 

In a previous paper (Valdarnini 2002) results from a set of hydrodynamical SPH sim- 
ulations of galaxy clusters have been used to investigate how final X-ray properties of the 
simulated clusters depend upon the simulation numerical resolution and the chosen star for- 
mation (SF) prescription. For a chosen SF model, final X-ray luminosities have been found 
to be numerically stable and consistent with data, with uncertainties of a factor ^ 2. In these 
simulations the metal enrichment of the ICM has been included with a minimal number of 
prescriptions and it was shown that for the simulated clusters, the final iron abundances of 
the ICM are in broad agreement with measured values. 

Chemical evolution in hydrodynamical SPH simulations has already been considered in 
a variety of contexts (Steinmetz & Miiller 1994; Raiteri, Villata Sz Navarro 1996; Carraro, 
Lia & Chiosi 1998; Buonomo et al. 2000; Mosconi et al. 2001; Churches, Nelson & Edmunds 
2001; Lia, Portinari & Carraro 2002; Aguirre et al. 2001). Metzler & Evrard (1994) have 
investigated the metal enrichment of the ICM in P3MSPH simulations of galaxy clusters 
in a standard cold dark matter (CDM) scenario. The authors have not included radiative 
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cooling in the simulations and have adopted a phenomenological prescription to model the 
chemical enrichment. 

In a more recent paper, Kravtsov & Yepes (2000) have estimated SN heating of the ICM 
using fixed-grid Eulerian hydrodynamical simulations. However, the numerical resolution of 
the simulations was not adequate to study the evolution of simulated clusters. The number 
of SNe occurring in a given cluster was estimated statistically from many small-box galaxy 
formation simulations. The implementation of a self-consistent metal enrichment model for 
the ICM in hydrodynamical simulations is therefore important for investigations of the ICM 
metal evolution. This is particularly relevant in connection with the observed metallicity 
gradients and for assessing the reliability of the calculated amount of ICM heating from 
SNe, inferred from measured metallicity abundances. 

The main aim of this paper is to analyze, in hydrodynamical SPH simulations of galaxy 
clusters, the dependence of the final iron abundance on a number of model parameters that 
control the ICM metal enrichment. This is done in order to obtain for the simulated clusters 
a final ICM distribution which can consistently fit a set of observational constraints, such 
as the observed iron abundances and at the same time the luminosity-temperature relation. 
Implications for the ICM heating from SNe are also discussed. This paper constitutes a gen- 
eralization of a previous work (Valdarnini 2002, hereafter V02), where the investigation was 
mainly concerned with the numerical stability of different SF models. Section 2 presents the 
hydrodynamical simulations of galaxy clusters that have been performed. Section 3 describes 
the star formation algorithm and the modeling of ICM metal enrichment implemented in 
the simulations. Results are presented in Section 4 and the conclusions in Section 5. 

2 SIMULATIONS 

Here I give a short description of the simulations. Further details can be found in V02. 
The cosmological model considered is a fiat CDM model, with a vacuum energy density 
= 0.7, matter density parameter 0,^ = 0.3 and Hubble constant h = 0.7 in units of 
100Kmsec~^Mpc~^ . The primeval spectral index of the power spectrum n is set to 1 and 
Qb = 0.015/i~^ is the value of the baryonic density. The power spectrum of the density 
fiuctuations has been normalized in order to match at the present epoch the measured 
cluster number density (Eke, Cole & Frenk 1996; Girardi et al. 1998). Initial conditions 
for the cluster simulations are constructed as follows. A collisionless cosmological N-body 
simulation is first run in a L = 200h^^ Mpc comoving box using a P3M code with 84^ 
particles, starting from an initial redshift Zin. At z = cluster of galaxies are located 
using a friend-of-friend algorithm, so to detect densities ~ 200Q^'^ times the background 
density within a radius r2oo- The corresponding mass M200 contained within this radius is 
defined as M200 = (47r/3)r2mPc^c'^200' where = 187f2~°-^^ for a fiat cosmology and pc is 
the critical density. The 40 most massive clusters within the simulation box are identified 
according to this procedure and the most massive and least massive cluster (labels 00 and 39, 
respectively) of this sample are selected for the hydrodynamical simulations. This procedure 
has been already followed in V02 and the initial conditions of the cosmological simulation are 
the same. In addition to these two clusters, two other clusters are considered. The original 
sample is enlarged to include 80 more clusters of decreasing virial mass. This amounts to a 
total of 120 clusters. The additional clusters extracted from the new sample are : the cluster 
highest in mass after L39 (L40) and the least massive of the sample (LI 19). Table 1 lists 
the properties of the four selected clusters. 

For each of the four test clusters, hydrodynamical TREESPH simulations are performed 
in physical coordinates. The initial conditions of the hydro simulations are determined as 
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Table 1. Reference values at z = for the four simulated clusters used in the numerical tests. M200: cluster 
mass within r200 in units of h~^MQ, r200 is in units of Mpc, ai is the central 1-D dark matter velocity 
dispersion in Kmsec~^ , Tm is the mass- weighted temperature in keV, fg is the ratio of the mass of gas 
within the radius r = 0.5h~^ Mpc to the total cluster mass within that radius. 
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Table 2. Numerical parameters of the simulations. gravitational softening parameter for the gas in 
Kpc. mass of the gas particles in h~^MQ. : mass of the dark particles, "^i number of gas particles 
inside the Lc/2 sphere at 2 = Zi^- : as in the previous column but for dark particles. ^: total number of 
simulation particles, including those in the external shell of radius Lc. ^ : initial redshift for the simulation 
If Zin > 20 the softenings are held fixed in comoving coordinates until z = 20, after which are kept fixed in 
physical coordinates. 

follows. The cluster particles at ^ = within r2oo of the test cluster are located in the original 

1 /3 

simulation box back at A cube of size Lc — 15 — 25 h~^Mpc oc Mjoq enclosing these 
particles is placed at the cluster centre. A lattice of grid points is set in the cube. At each 
grid point is associated a dark matter particle and a gas particle, of corresponding mass 
and coordinates. The particle positions are then perturbed, using the same initial conditions 
of the cosmological simulations. Finally, the particles for which the perturbed positions lie 
inside a sphere of radius Lc/2 from the cube centre are kept for the hydro simulations. 

To model the effects of the external gravitational fields, the inner sphere is surrounded 
out to a radius Lc by dark matter particles with a mass eight times larger than the sum of 
the masses of a gas particle and a dark matter particle of the inner sphere. For each particle, 

1/3 

gravitational softening parameters are set according to the scaling £j oc m/ . The numerical 
parameters of the hydrodynamical simulations are given in Table 2. 

The simulations with index LOO and L39 have a number of gas particles Ng ~ 22, 600. 
For this mass resolution, the corresponding runs in V02 have been found to give fairly stable 
final X-ray luminosities for the simulated clusters. Additional runs have been considered, 
with an increased resolution with respect to the standard runs. The high resolution (H) runs 
have Ng ~ 70, 000 and the very-high resolution (VH) runs have Ng ~ 210, 000. The other 
numerical parameters have their values scaled accordingly. 

The gravitational forces of the hydrodynamical simulations are computed using a hier- 
archical tree method with a tolerance parameter 6 = 1 and taking into account quadrupole 
corrections. The hydrodynamical variables of the gas are followed in time according to the 
SPH Lagrangian method (Hernquist & Katz 1989 , and references cited therein). In SPH, 
local fluid variables are estimated from the particle distribution by smoothing over a num- 
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ber of neighbors. A common choice is the i?2-sphne smoothing kernel Ws{r, h) (Monaghan & 
Lattanzio 1985), which has compact support and is zero for interparticle distances \r\ > 2h. 
The smoothing kernel is normalized according to J W{f, h)df = 1. The smoothing length 
h fixes the spatial resolution of the simulation. The resolution is greatly increased in high 
density regions when individual smoothing lengths hi are allowed, so that the number of 
neighbors of a gas particle is nearly constant 32). A lower limit to the smoothing lengths 
hi is set by hi > hmin = ^c//4, where Eg is the gravitational softening parameter of the 
gas particles. The time integration is done allowing each particle its own time step. The 
accuracy of the time integration is controlled by a number of constraints that the individual 
time steps must satisfy ( see, e.g., Valdarnini, Ghizzardi & Bonometto 1999). The minimum 
allowed time step for the gas particles is 6.9 ■ lO^yr. 

The thermal energy equation for the gas particles includes a term which models the 
radiative processes of an optically thin plasma in ionization equilibrium. The total cooling 
function Ac depends on the gas temperature and metallicity. The cooling function takes into 
account contributions from recombination and collisional excitation, bremsstrahlung and 
inverse Compton cooling. Heating from an ionizing UV background has not been considered. 
The cooling rate of the gas in the simulations is then dependent on the gas metallicity. This 
is an important difference with respect the previous simulations (V02) and is essential in 
order to analyze low-temperature {Tx — 2keV) clusters consistently. The dependence on 
the metallicity has indeed larger effects, as it increases the cooling rate. The small scales 
resolved by the simulations will cool faster than larger ones. The inclusion of cooling with 
its metallicity dependence has then a strong impact also on the formation of larger clusters 
because of the hierarchical growth of structure. 

Tables of the cooling rates as a function of the temperature and gas metallicities have 
been constructed from Sutherland & Dopita (1993) and stored in a file. During the simu- 
lations a cubic spline interpolation is then used to calculate from the tabulated values the 
cooling function Ac(T, Z) of a gas particle of given temperature T and metallicity Z. Here Z 
is the mass fraction of metals of the gas particle. Conversion from the metallicity Z to the 
corresponding value of the iron-to hydrogen ratio [Fe/H] {[X] = IoqiqX — IoqioXq) is done 
as in Sutherland & Dopita. A good fit to this relation is given by Eq. 1 of Tantalo, Chiosi 
& Bressan (1998). The X-ray luminosities are computed from the gas emissivities within 
the cluster virial radius according to the standard SPH estimator (see Eq. 8 of V02). The 
X-ray emissivity e{r) associated with a gas particle is calculated with a Raymond-Smith 
code (1977) as a function of the gas temperature and metallicity. 

3 STAR FORMATION AND ICM ENRICHMENT 

Cold gas in high density regions will be thermally unstable and subject to SF. In SPH sim- 
ulations SF processes have been implemented using a variety of algorithms. Here conversion 
of cold gas particles into stars is performed according to Navarro & White (1993). In V02 
it has been found that for this SF method final profiles of the simulated clusters are robust 
against the numerical resolution of the simulation. According to Navarro & White any gas 
particle in a convergent flow and for which the gas density exceed a threshold. 



will be in a collapsing region with its cooling time smaller than the dynamical time and 
is eligible to form a star particle. If these conditions are satisfied, SF will occur with a 
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characteristic dynamical time scale Td = y3'iT/lQGpi. The probability that a gas particle 
will form a star in a time step At is then given by 

p=l-exp{-At/Td). (2) 

A Monte Carlo method is used at each time step to identify those gas particles which form 
star particles. For these gas particles, a star particle is created with half-mass, position, 
velocity and metallicity of the parent gas particle. The star particle is decoupled from the 
parent gas particle when it is created and is treated as a collisionless particle. A gas particle 
is converted entirely into a star particle when its mass falls below 5% of the original value. If 
the cooling function Ac depends also on the gas metallicity Z, the density threshold criterion 
must be modified in order to take into account the increased cooling rate when Z > 0. The 
condition that in Navarro & White (1993) defines the gas density threshold pc,g now reads: 

_ 3{pc,g/fimp)kBT _ 
Tcool- ^^{T,Z) 

where TcooI is the cooling time, T ~ 10^'^ °K, fi = 0.6, ks is the Boltzmann constant and rrip 
is the proton mass. Because of the metallicity dependence of the cooling function Eq. (3) 
must be solved numerically. A plot of pc,g as a function of Z shows that pc^g{Z) is nearly flat 
for Z 10"'^ and has a fast decay above Z ~ 5 ■ 10"'^. 

It has been found that a good fit to Pc,g{Z) is the following approximation: 

p,,g{Z) ^ 7 ■ 10-2Vcm-V(l + 0.3(Z/Z_l,5))=^/^ (4) 

where 2'_i.5 ~ 1.28 ■ 10~^ is the metallicity corresponding to [Fe/H] = —1.5. 

Once a star particle i is created at the time ts{i) it will release energy into the surrounding 
gas through SN explosions. SN of type II (SNII) originate from the explosions of stars of 
mass rriu > m > 8Mq at the end of their lifetime T{m), here r(m) is defined as in Navarro 
and White (1993). 

Each SN explosion produces Esnii = £sn — 10^^ erg , which is added to the thermal 
energy of the gas, and leaves a ^ IAMq remnant. The number of SNII explosions associated 
with the star particle i in the time interval [t — At, t] is determined as 

ANii{t) = nisii) (j){m)dm, (5) 

where 0(m) is the IMF of the stellar population, is the root of T{m), t* = t — ts{i) = 
T, mgii) is the mass of the star particle and rriu is the upper limit of the IMF The 
normalization of the IMF is set to /q"|" (pmdm = 1, with this normalization Np^p = ms{i) 
is the number of stellar populations of the star particle i. Several forms of the IMF have 
been considered. A Miller-Scalo (1979) has been chosen for a consistent comparison with 
the simulations of V02. For this IMF rriu = 100. A standard IMF is the one of Salpeter 
(1955), where </)(m) oc m~^^~^^\ with x = 1.35. Finally, a less steep IMF is given by Arimoto 
& Yoshii (1987) for elliptical galaxies, for which a; ~ 1. The latter two IMF have rriu = 40. 
The index of the simulations with different IMF are presented in Table 3. 

The energy produced in the time interval [t — At, t] by the SN explosions of type II 
associated with the star particle i is EsnANh. This feedback energy is returned entirely to 
the nearest neighbor gas particles of the star particle i. The velocity field of the neighboring 
gas particles is left unperturbed by the SN explosion, since the typical SPH simulation 
resolutions are much larger than the size of the shell expansion (Carraro, Lia & Chiosi 

* Hereafter masses are in solar units 
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Table 3. Model parameters of the simulations. shape of the IMF; MS=Miller-Scalo,S=Salpeter and 
A=Arimoto-Yoshi. the cooling function includes also the dependence on the gas metallicity. : SN of 
type la are considered as additional sources of the ICM enrichment, a non zero entry is the value of the 
normalization constant (Eq. 7). For these simulation Wz = B2 and esN = 1 (see Table 4). 



1998). The energy is smoothed among the internal energies of the gas particles according to 
the SPH smoothing prescription. The internal energy increment Auj of the gas particle j is 
then 

where rrij is the mass of the gas particle, pj is the gas density, Wg is the B2 SPH smoothing 
kernel, 2/z| is the radius of a sphere surrounding A^^ ~ 32 gas neighbors of the star particle 
i and Vc = J2j Ws{rj — r^, hl)mj/ pj is a normalizing factor. Vc has been introduced to avoid 
that Aujrrij 7^ Esn^Nu. The smoothing length is constrained by the upper limit 

< SOKpc = Hm and there is no lower limit, unless explicitly stated. In high density 
regions the SN energy is then returned to the gas neighbors without being affected by the 
SPH resolution length. The choice of the value of Hm is discussed in sect. 4.2, where the 
dependence of the simulated cluster profiles on a number of parameters is analyzed. 

The number of SN of type la (SNIa) associated with the star particle i in the time 
interval [t — At, t] has been determined according to Greggio & Renzini (1983). The scheme 
implemented follows Lia, Portinari & Carraro (2002). The number of SNIa events at the 
epoch t is 



NsNiai< t) = m,{i) A 



f^b.sup 


' 1-0.5 


/ 0("^f)) 


/ fifJ')dp 




J fJ-inf 



drUb, (7) 



here rrih is the mass of the binary system, < p = 7722/7725 < 0.5 is the mass fraction 
of the secondary and f{p) = 24/i^. The mass rrih of the binary system lies in the range 
3 = mh inf < nil, < ^b,sup = 12 = 2M„p, 7772 = 77i2(r) is the mass of the secondary which 
ends its lifetime at the time t = T{m2) + tgii) and p^j = M AX {m2{T) / nih, 1 — Mup/nii,}. 
The lower limit to the mass of the secondary, 0.9 < 7712, is given by the age of the universe. 
The normalization constant A is set to A = 0.07 in order to match the estimated rate of type 
la SNe in galaxies (Portinari, Chiosi & Bressan 1998). Inverting the order of integrations 
one obtains: 



NsNiai< t) = ms{i) A / 24m: 

•^"12 (r) 



2 



3 dnib 



dm2, (8) 



where rrib^rnin = M AX {mb,inf , "21122} , mb^max = MIN{mb,sup,rn2 + M^p) for the number 
of SNIa explosions associated with the star particle i at the age t. The number of SNIa 
explosions in the time interval [t — At,t] is then ANia = {NsNia{< t) — NsNia{< t — At)). 
The corresponding explosion energy is esNia^Nja, with esNia — ^sn (Woosley & Weaver 
1986). This energy is smoothed as in Eq. (6) over the nearest neighbor gas particles of the 
star particle i. 
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run 


IMF 


SNIa 


h " 
A/ 




e 

SN 


h"^ ■ 

rmn 


A 


A 


0.07 


30 


B2 


1 





A2 


A 


0.07 


10 


const 


1 





A3 


A 


0.07 


30 


const 


1 





A4 


A 


0.07 


30 


const 


1/10 





A5 


A 


0.07 


30 


const 


1 


12 


A6 


A 


0.07 


30 


const 


1 


6 



Table 4. Additional model parameters for runs with AY IMF (label A). maximum smoothing length for 
the metal enrichment spline. Functional form of the chosen spline for the metal enrichment of the gas; 
B2'- standard SPH, const for a uniform spread. : SN explosion energy in units of 10^^ erg. : minimum 
smoothing length for the metal enrichment. 



SN explosions also inject enriched material into the ICM, thus increasing its metallicity 
with time. The mass of the k — th heavy element produced in a SN explosion is defined 
as the stellar yield ys,k{'m)i with S'=II or la. The total yield is the sum over the masses of 
the heavy elements : ysijn) = Hkys,k{jn)- For SNIa the yield yia,k is a constant which is 
independent on the progenitor mass. The adopted SNIa yields are those of Iwamoto et al. 
(1999, W7 model), with yia = 1-4 and yia,Fe = 0.7. The yields yij,k of type II SN depend on 
the progenitor mass m and it is useful to define an average yield 

< yii,k >= ,/ • 9 

The upper limit m„ is defined according to the IMF and m/ is the lower bound for SNII 
progenitors. Here it is assumed mi = 8. For type II SNe the predicted theoretical yields suffer 
from a number of uncertainties. According to the assumed massive star physics the yields 
of different models can differ by a factor ~ 2 or more. A detailed discussion can be found in 
Gibson, Loewenstein & Mushotzky (1997). The theoretical model chosen here is model B of 
Woosley & Weaver (1995). For this model the explosion energy of massive stars is boosted 
and the produced metals are not reabsorbed in the SN explosion. The yields of this model are 
therefore the most favorable from the point of view of the amount of iron synthesized from 
SNII. Final metal abundances for a different theoretical SN explosion model can be obtained 
from the SNII component of the gas metallicity profile by a rescaling of the yields. This is 
a valid approximation as long as there are small differences between the yields yu which 
determine the gas metallicity Z, and therefore the local cooling rate and the SF threshold 
(Eq. 4). 

Woosley & Weaver (1995) yields are available down to = 11. The yields between 
mi = 8 and m = 11 have been obtained with a linear extrapolation, the associated errors 
are however small since SN nucleosynthesis is negligible below m = 11. From Table 1 of 
Gibson, Loewenstein & Mushotzky (1997) < yn,Fe >— 0.14 for the Woosley & Weaver B 
model with a Salpeter IMF, m/ = 10 and m„ = 50. I obtain the same value for an Arimoto- 
Yoshi IMF with m/ = 8 and m„ = 40. The mass of heavy elements which has been injected 
into the ICM at the age t due to the SNII explosions associated with the star particle i is 
determined as (Tinsley 1980) : 



'ms,zii{< t) = ms{i) < [m-mr - yii{m)] Zs{i)(j){m)dm + / yji{m)(f){m)dm > , (10) 

\Jm{T) Jm{T) \ 

where Zs{i) is the metallicity of the star particle i , m(r) is the SNII progenitor mass at the 
epoch T = t — ts{i) and m^ ^ 1.4 the remnant mass. Similarly, the total mass returned to 
the ICM is 
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m 



(< t) = ms{i) / [m — rrir] (p{m)dm. 



J m(T) 



(11) 



The ejected masses from star particle i in the time interval \t — At, t] are then Arris^esp = 
fTT's,esp{t) — ^s,esp{t — At), with esj) = ZII or ejec. For each star particle i, these masses are 
calculated at each time step according to the stellar particle age r. The Arris^esp are then 
distributed over the gas neighbors of the star particles, using the same smoothing procedure 
adopted (Eq. 6) for returning the SN feedback energy to the gas, i.e. 



where Amzii{j) is the mass increment of the metallicity of type II of the gas particle j, 
neighbor of the star particle i, and Wz is the smoothing kernel. The same relation holds for 
the total ejected mass Am^^ejec- This smoothing procedure is similarly applied to calculate 
the metallicity increment of the gas particles due to SNIa ejecta. Since yia is a constant for 
SNIa the mass increment, Arris^ejec is defined as Arris^ejec = ANjait) yia- Each gas particle, 
as well as a star particle, has two different metallicity variables. These variables represent, 
respectively, the total mass of metals originated from SNII and SNIa events. This is necessary 
in order to discriminate the different contributions of SNII ejecta from that of SNIa in the 
final profiles of the gas metallicity. The above integrals are calculated numerically as a 
function of the progenitor mass m = T~^{t) and stored in tables. A linear interpolation 
procedure is used during the simulation to obtain from the tables the values of the integrals 
at the age t. 

Since the choice of the smoothing kernel Wz is somewhat arbitrary, a standard choice is 
the SPH smoothing kernel (Mosconi et al. 2001; Lia, Portinari & Carraro 2002). However, if 
the ejected material is deposited at the end of the shock expansion, a less steep deposition 
profile is a better description of the deposition mechanism. For a uniform distribution, 
Wz oc const. Aguirre et al. (2001) have investigated the metal enrichment of the diffuse 
intergalactic medium in cosmological simulations. They have considered a power-law shape 
W{r, h) oc r" for the deposition kernel. They assumed a default value a = 3, the uniform 
case considered here corresponds to their a = 2. The simulations have been performed using 
for Wz the SPH smoothing kernel as the default kernel. Table 4 reports the index of the 
runs for which Wz = const. Each simulation has a label obtained by merging together the 
labels of Table 2 and 3, or Table 2 and 4, which determine the simulation parameters. 



A measure of the regularity of the numerical cluster sample can be assessed from a com- 
parison of the mean gas fraction of the four test clusters fg = Mg/M^ within the radius 
r = 0.5h~^Mpc, with the corresponding values of the Matsumoto et al. (2000) sample (Fig. 
1). This sample has been chosen for a comparison of the global cluster iron abundances with 
the ones predicted by several runs (see below). The sample values of fg shown in Fig. 1 have 
a large scatter around the cosmological value of the model = flb/Qrn = 0.1, while those of 
the simulated clusters are systematically smaller. This indicates that the simulated clusters 
are regular objects (Evrard, Metzler 8z Navarro 1996) in a quiet dynamical state, without 
having undergone recent mergers, a conclusion supported also from previous substructure 
analyses of X-ray maps (Valdarnini, Ghizzardi & Bonometto 1999) for a ACDM cluster 
sample with the same cosmological initial conditions. Therefore, the profiles of simulations 
with different model parameters can be consistently compared without biases that can follow 
from the presence of substructure. 




m,- 1 



(12) 



4 RESULTS 
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Figure 1. Baryonic fraction for tlie Mastumoto et al. (2000) sample of nearby (z < 0.1) clusters, fg = Mg/Mx is evaluated 
at r = Q.bh~^ Mpc. Open symbols refer to the simulated clusters. The dashed line is the cosmological value (Q^/Qm) for the 
assumed model. 

4.1 Simulations with different IMF and cooling 

Results from simulations with different IMF and cooling parameters are discussed first. 
The relevant parameters for these runs are presented in Table 3. The simulation results are 
shown here only for the cluster ACDM39; for the other three test clusters there are not 
qualitatively relevant differences in the final profiles. The simulations have been performed 
keeping the numerical resolution, which is given by the index L39 of Table 2, fixed. The 
first three runs of Table 3 (L39NZlal, L39NZla2 and L39NZla3) correspond to the following 
choices of the IMF : Miller-Scalo (MS), Salpeter (S) and Arimoto-Yoshi (AY), respectively. 
All the other parameters have been left unchanged. These simulations do not consider the 
possible dependence of the cooling function on the gas metallicity and the contribution of 
SN of type la to the gas metal enrichment. For this choice of parameters, the run L39NZlal 
is just case cl39-10 of V02, in order to compare final profiles with previous findings. L391a4 
has the metallicity dependence of the cooling function switched on. The other two indexes 
of Table 3 (S and A) are for runs with SNla as a heating source of the ICM, as well as of 
metal enrichment. 

Fig. 2 shows the final profiles of temperature, gas and iron abundances in solar units 
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Figure 2. For the cluster ACDM39, a comparison of simulations with different choices of the coohng function and the IMF. 
Bottom left: SF rates as a function of time. The other panels refer to z = 0. From the top left : temperature, gas densities and 
iron abundances profiles versus radial distance. Iron abundances arc in solar units. In the density panel the lower profiles are 
those of the metals. The meaning of the simulation labels is given in Table 2 and 3. 



for the different runs. The SF rate (SFR) as a function of age is also plotted. While a 
comparison of the iron abundances obtained from simulations with available data will be 
discussed later, several conclusions about the choice of the IMF can already be drawn from 
the iron profiles seen in Fig. 2. For a cluster like ACDM39, the iron abundance expected at 
a distance ^ OAMpc (~ 0.2r2oo) from the centre is in the range ~ 0.3 — 0.35 (De Grandi & 
Molendi 2001). Hereafter iron abundances are given is solar units {Fe/H = 4.68 ■ 10"^ by 
number). The iron profiles of Fig. 2 show that a Miller-Scalo IMF is completely ruled out as 
a possible IMF and cannot produce the amount of iron required by observations. The same 
conclusion holds also for a Salpeter IMF, for which Fe/H ^ 0.1 at r ~ OAMpc. The only 
IMF for which the simulation gives a significant amount of iron is Arimoto-Yoshi (AY). The 
iron profile of the S run L39NZIa2 is always a factor two below that of L39NZIa3 (AY) . At 
r ~ OAMpc the latter profile gives Fe/H ~ 0.2 — 0.15. These values are still below the ones 
required by observations, nevertheless it appears that the AY IMF gives the best results 
from the point of view of the amount of iron required to fit the data. These conclusions are 
in agreement with those of Loewenstein & Mushotzky (1996). 
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These iron profiles are originated from tlie metal ejecta of SN II in the ICM, with the 
bulk of the SN explosions that has enriched the ICM already at 2; ^ 0.7. Therefore, the 
gradients in the final iron profiles indicate that the cluster was already dynamically relaxed 
at this epoch, without major mergers of substructure which could have remixed the ICM 
and erased the original gradients. The most important differences in the profiles arise when 
the cooling function also depends on the gas metallicity (L39NIa4). For this run, the SFR is 
much higher than in the previous cases. For the range of temperatures considered here, this 
is a direct consequence of the higher cooling rate. As a consequence, the density threshold 
criterion for SF is lower and the SF activity is much higher than in the no-metal cases. The 
top right panel of Fig. 2 shows that final gas density profiles are not strongly affected by the 
metallicity dependence of the cooling rate. 

One of the most important consequences of considering metallicity effects for the cooling 
rate is seen in the temperature profiles. From the top left panel of Fig. 2 there is clear evidence 
that the temperature profile of L39NIa4 has a different shape from the other profiles already 
at r ~ 2{]{]Kpc from the cluster centre. Between this distance and the cluster centre, the 
profile shows values of the temperature higher than in the runs with the no-metal cooling. 
The metallicity dependence of the cooling rate affects the final temperature profiles at the 
cluster centre through two main effects. The first is the term Ac(T, Z), which enters in the 
gas thermal equation. Even for gas temperatures above ~ 2keV the increase in the cooling 
rate with respect the no-metal case can be significant for high metallicities {Fe/H ^ 1), 
which can be present at the cluster centre if there are strong metallicity gradients. The 
second effect is indirect; because of the increased cooling rate, SF has been higher in the 
past than in the no-metal case, and as a consequence there is a larger amount of cold gas 
that has been converted into stars and a deeper potential well at the cluster centre. As a 
result , there is a higher inflow at the cluster centre of the surrounding high-entropy gas than 
in the runs that do not consider the metallicity dependence of the cooling function. This in 
turn implies higher temperatures at the cluster centre (see sect. 4.2). An important result 
that therefore follows from these simulations is that the temperature profiles in the cluster 
inner regions cannot be considered as approximately flat. This is in disagreement with what 
has been found in V02, where the simulations included radiative cooling and SF, but did not 
take into account the metal contribution to the cooling. These conclusions are rather general 
and are valid not only for ACDM39, which has a mass- weighted temperature Tm — 3keV, 
but also for a cluster like ACDMOO, for which Tm — 6keV (see Fig. 5 below). 

The iron abundance profile of the run L39NIa4 is lower by a factor ^ 50% with respect 
that of L39NZIa3, where the cooling function does not depend on the gas metallicity. This 
is also confirmed by contrasting the metal density profiles in the top panel of Fig. 2. If the 
cooling rate increases when the metallicity dependence is taken into account, then the final 
amounts of stars and metals at the cluster centre are expected to be higher than in the 
runs with the no-metal cooling. This is in contrast with what is found. It is not clear how 
this anti-correlation between the final metal ICM abundance and the metal dependence of 
the cooling rate originates. A possible explanation lies in the strong SF activity that occurs 
earlier for the run under consideration. As a consequence, most of the ICM is metal enriched 
by stars at earlier epochs. An iron profile that is in better agreement with data is recovered 
when SNe of type la (L39A) are also considered as sources of metal enrichment for the ICM. 
For this run, the iron abundance at ~ SOOKpc is ~ 0.2 , about two times that of the parent 
run with the same parameters but without SNIa (L39NIa4). This is obtained for A = 0.07 
(Eq. 7). According to Portinari, Chiosi & Bressan (1998) this choice of A corresponds to a 
number of SNIa ^ 20% of the total number of SNe in the Galaxy. A numerical simulation 
with A = 0.04 has shown that the amount of iron that is injected from SNIa in the ICM 
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Figure 3. For the clusters ACDMOO and ACDM40 final iron profiles from SN Il/Ia are displayed for simulations with a Salpeter 
IMF (see Table 3). Filled symbols are for the corresponding very-high resolution runs. 

is roughly proportional to A. The above results indicate that for an AY IMF as much as 
^ 50% of the iron in the ICM comes from SNIa, provided that the current rate is normalized 
to match that estimated from spiral galaxies. 

The stability of the final iron profiles against the simulation's numerical resolution can 
be estimated from Fig. 3, where for a chosen IMF (S) the iron profiles of the two clusters 
ACDMOO and ACDM40 are displayed for the standard resolution and the very-high resolu- 
tion runs (VH). There is a good agreement between the profiles with a different resolution, 
with small differences localized within r ^ O.lMpc. 

4.2 Simulations with different metal ejection parameters 

Values of the global iron abundances Ap^ = Mpe/Mu are listed in solar units in Table 5 for 
the various simulations. These values have been calculated at a fiducial radius r = 0.5/i~^Mpc 
and are in the range ~ 0.15 — 0.2 for the simulations with index A. From the Matsumoto et 
al. (2000) sample of nearby clusters [z < 0.1) the estimated abundances at 0.5h~^Mpc are 
all above ~ 0.2 and in the range ^ 0.2 — 0.4 (see later Fig. 5). Therefore, it is important to 
investigate the effects on the final iron abundances of adopting prescription of metal ejection 
different from that of the A runs, keeping fixed the other parameters of the simulations. 

The largest amount of iron ejected in the ICM are obtained for the AY IMF, and this is 
the IMF that is chosen in the remaining of the paper in order to account for the measured 
iron abundances. Hence for this IMF, different model parameters of metal ejection have been 
tested. Table 4 presents the index of the runs and the associated relevant parameters. 

The metal enrichment of the gas is modeled according to Eq. 12, with the mass that 
is ejected by a star particle i in the interval At distributed over the gas neighbors, the 
mass fraction is weighted according to the smoothing kernel Wz{r,h). For the A runs the 
kernel Wz is the standard SPH B2 spline. The radial profile of the deposition kernel is 
one of the parameters that governs the distribution of the ejected metals among the gas 
particles surrounding the star particle. According to Aguirre et al. (2001) the form of the 
distribution function Wz{r,h) can be generically assumed with a radial power law oc 
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Figure 4. For the four test clusters final iron profiles arc plotted for runs with different model parameters of metal enriehment 
(Table 4) and numerical resolution. 

behavior. A radial profile shallower than that of the B2 spline clearly implies that more 
metals are deposited near the limiting radius. Therefore, gas particles that are not part of a 
SF activity can be metal enriched and diffuse metals into the ICM. The metal enrichment 
of the ICM is expected to be higher in this case because metals are more effectively mixed. 
This has been also pointed out by Mosconi et al. (2001). Alternatively to the B2 spline, the 
shape of the deposition profile chosen here is Wz = const, which corresponds to the case of 
a uniform metal distribution. 

The simulations with index A3 in Table 4 have Wz = const and are the mirror simulations 
of the runs with index A, the only different parameter being the choice of the deposition 
kernel. Final iron profiles are compared against those of models A in Fig. 4. The four panels 
of Fig. 4 refer to the four test clusters. In each panel final iron profiles are plotted for the 
runs with different metal ejection parameters and different numerical resolutions. The plots 
of Fig. 4 show that, for models A3, the profiles are higher than that of models A by a factor 
that is in the range ^1 — 2. The differences are largest for ACDMOO and ACDM39, modest 
for ACDM40. The global iron abundances Ap^. are now in the range ~ 0.23 — 0.35 (Table 
5), in better agreement with data than the values of models A. These results demonstrate 
that the choice of the deposition profile is a key parameter in determining iron abundances 
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in the ICM. The differences in the final profiles between the runs A3 and A gives a measure 
of the scatter in the final abundances that is inherent to the choice of different kernels. 

It must be stressed that other choices are clearly possible; in particular a steeper IMF may 
even require a deposition kernel with a positive radial derivative. A quantitative comparison 
with several data of the profiles obtained from model A3 for the four test clusters is performed 
later. The results show that for models A3 the simulations are in good agreement with the 
measured values and the parameters of the A3 runs are then taken here as those of a 'fiducial 
model' against which to compare the results from the other simulations. 

The stability of the simulated iron profiles against the numerical resolution is tested 
in Fig. 4 by plotting (filled symbols) profiles of high-resolution simulations (H in Table 2) 
against the corresponding ones with standard resolution. As can be seen, there are no large 
differences in the final profiles between the high and standard resolution runs. The only 
important exception is for the coldest cluster ACDM119 (T^ ~ LbkeV). The iron profile 
of L119HA3 is much lower than that of the standard resolution run L119A3. This strongly 
suggests that there are numerical resolution problems that affect the profile of ACDM119, 
when the numerical resolution is that of the standard runs. Increasing the numerical resolu- 
tion implies that smoothing lengths are smaller. As a consequence the metals are distributed 
among gas neighbors in smaller volumes, this implies a smaller mixing of metals in the ICM 
(see also Mosconi et al. 2001). Clearly, this effect should not be seen (cf. Fig. 3) and if present 
implies that low-resolution runs are undersampling the gas distribution. In order to reliably 
predict metallicity profiles it is then safe to assume that simulations of low-temperature 
(T ^ 2)keV) clusters require at least the mass resolution of high-resolution runs. To explain 
the low values of the Ape for the run L119HA3 an alternative hypothesis is that the metal 
enrichment of the ICM by SNe is characterized by a minimum diffusion length. Simulations 
have been performed without a lower limit h^^^^ for the smoothing lengths /i| of the diffu- 
sion kernel Wz of ejected metals. For cool clusters high resolution runs may therefore imply 
values of hi in the central high density regions below /i^j^. This in turn would imply lower 
values for the metallicity profiles because of the reduced mixing of the ejected metals in the 
ICM. This effect is very similar to the behavior expected when the numerical resolution is 
insufficient and high-resolution runs yield lower profiles than those of low-resolution runs. 
The profiles of Fig. 3 show that this effect is relevant only for the less massive, cool, clusters 
( ^ ?)keV). For the numerical parameters of the high-resolution runs (H) this implies that 
must be of the order of ~ few Kpc. In order to investigate the dependence of the final 
metallicity profiles on the value of /i^in simulations have been performed for the cluster 
ACDM119. Models A5 and A6 have the same parameter of model A3, but with = 12 
and 6Kpc, respectively. Iron profiles for these two runs are plotted in the panel of Fig. 4 for 
the cluster ACDM119. The results demonstrate that requiring a minimum diffusion length 
is highly effective to obtain for cool clusters a large-scale mixing of iron in the ICM. The iron 
abundances Ap^ are now in the range ~ 0.4 — 0.5, in better agreement with the observational 
values for cool clusters. 

Another parameter which is important to determine the amount of metal enrichment 
of the ICM is the upper limit Hm (see Table 4) of the smoothing lengths h^. These are 
defined in Eq. 6 as the half-radius of a sphere enclosing ~ 32 gas neighbors of the generic 
star particle i. The standard value assumed for the simulations is = >iOKpc ^ . The 
sensitivity of the final profiles to the assumed values can be estimated from the plots of 
Fig. 4. For the runs A2, Hm = lOKpc. Contrasting the profiles with those of the A3 runs 

t This is the same value used in V02; the value quoted in the text of V02 is erroneously half the true value 
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runs Ap^{0.5) 



ACDMOO LOOS/LOO VHS/L00A/L00HA/L00A2/L00A3/L00HA3 0.106/0.145/0.236/0.251/0.147/0.365/0.376 

ACDM39 L39A/L39A2/L39A3/L39HA3 0.188/0.127/0.288/0.229 

ACDM40 L40S/L40VHS/L40A/L40HA/L40A2/L40A3/L40HA3 0.08/0.106/0.208/0.203/0.185/0.277/0.245 

ACDM119 L119A/L119A3/L119HA3/L119A4/L119HA4/L119HA5/1119HA6 0.186/0.377/0.181/0.392/0.147/0.475/0.283 

Table 5. Average value in solar units (4.68 ■ 10~^) of the iron abundance Aj?^ = Mpg{< r)/M[j{< r) for 
the considered runs. Ap,, is evaluated at r = 0.5h~^ Mpc. 

shows clearly that models A2 yield final abundances well below those of models A3 and 
comparable to the abundances of models A. For the runs A2, the abundances Ape are of 
the order of ^ 0.1 — 0.13, a factor ^ two smaller than the values of runs A3. The value of 
Hm = lOKpc is clearly ruled out by the measured iron abundances. The choice Hm = 30Kpc 
gives much better results for the metal distribution in the ICM. In fact, the values of hf are 
rarely fixed by this upper limit. A radial binning of the squared distribution {hlY shows 
that on average the rms values of the smoothing lengths are below lOKpc in high-density 
regions (r O.lMpc), and grow up to ~ 20Kpc in low-density regions (r ~ 500Kpc) for 
which p/pc — 10. This shows that the choice of the value of hM has a relevant impact on 
the final metallicity profiles at radial distances r ^ O.lMpc, in low-density regions where is 
most likely that gas particles that are not in a SF region can get metal enriched and diffuse 
metals. These values of the smoothing lengths correspond to maximum radii ~ AOKpc and 
are much higher than the upper limit of ~ lOKpc estimated by Ezawa et al. (1997) for the 
diffusion of ions in a Hubble time in a low-density plasma [p ~ lOpc) at a temperature of 
~ 4keV. On the other hand the Ezawa et al. limit corresponds here to hM = 5Kpc, as it 
has been found for models A2 this choice of hM would imply for the simulated clusters final 
iron abundances in the ICM much lower than the measured values. The plots of Fig. 3 & 4 
show also that for clusters with temperatures above ^ 2keV the profiles of model A3, with 
hM = SOKpc, do not depend on the numerical resolution of the simulations. A possible way 
of reconciling these discrepancies lies in the fact that analytical estimates of the maximum 
diffusion length of metals in the ICM do not take into account the mixing of metals that can 
occur because of the dynamical interactions between cluster galaxies and the ICM (Dupke & 
White 2000b). Simulation profiles show that to model the mixing of metals the best results 
are obtained for the parameters of model A3. 

4.3 Comparison with data 

For the four test clusters, observational variables from simulations with different model 
parameters are compared in Fig. 5 against a number of data. For model A3 the projected 
emission weighted temperature profile are plotted as a function of radius in panel (a). Data 
points are the mean error- weighted temperature profiles of 11 cooling fiow clusters from 
De Grandi & Molendi (2002). The profiles have been calculated according to Eq. A3 of De 
Grandi & Molendi. The cluster centre is defined as the maximum of the gas density, and a 
peak emission criterion for defining the centre does not modify the calculated profiles in a 
significant way. Each smoothed profile has been rescaled to match the last data point. There 
is a remarkable agreement of the simulated profiles with data, the only important exception 
being for the two innermost bins. The temperature profiles show a radial increase toward 
the cluster central region, followed by a strong drop at the cluster centre. This feature is 
common to all the runs and is not shared by the data, for which the innermost bin has a lower 
temperature than the nearby ones. This behavior of the temperature profile is robust and 
is not sensitive to an increase of the numerical resolution. On the other hand, this increase 
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of the cluster temperature at the centre is a consequence of the entropy conservation during 
the galaxy formation and the subsequent removal of low-entropy gas (Wu & Xue 2002a). A 
possible explanation for this discrepancy lies in the fact that what is being measured are 
spectral temperatures. Mathiesen & Evrard (2001) argue that spectral fit temperatures are 
weighted by the number of photons, so line cooling from small clumps can bias the spectrum 
toward lower temperatures. This issue can be settled performing a spectral fit analysis as in 
Mathiesen & Evrard (2001), whose simulations did not include radiative cooling. 

The projected metallicity profiles are displayed as a function of distance in panel (b). 
Data points represents the mean profile from the nine cooling flow clusters of De Grandi 
& Molendi (2001). For model A3, the iron profile of L00A3 is the one in better agreement 
with data. The profile of L39A3 also has a shape similar to that of L00A3, but with a much 
higher abundance at r ~ 0.02r2oo ( Fe/H ~ 0.7 — 0.8). The profile of L40A3 agrees fairly 
well with data only for the first three radial bins. The overall shape is similar to that of the 
other two runs, but with a lower amplitude. The systematic difference between the profiles 
of runs L39A3 and L40A3 is mostly due to the different dynamical histories of the two 
clusters. Therefore, there are final uncertainties in the final profiles which can be as high 
as ~ 50% and are related to the cluster dynamical evolution. At outer radii all the profiles 
show a radial decay steeper than that of data points, which in fact can be considered to have 
an almost constant profile. It appears very difficult to modify the model parameters of the 
simulations in a way such that the simulated profiles match the data points at outer radii, 
without also increasing central abundances. 

The flatness of the observed profiles suggest that at, early epochs, a series of merger 
events has erased the existing metallicity gradients. Accordingly, the metal abundances from 
SNII do not show significant spatial gradients, while the iron abundance gradient can be 
attributed to SNIa (Dupke & White 2000b). In this scenarios the metal excess distribution 
is an indicator of the optical light distribution of early type galaxy, as it has been at least 
partially confirmed by De Grandi & Molendi (2001). This is not the behavior of the simulated 
cluster sample. Fig. 3 shows for two clusters the iron radial distribution originated from 
both SNII and la. The iron distributions are very similar, because the timescales of metal 
ejection is much higher for SNIa than for SNII this is indicative that the dynamical evolution 
of the two clusters has been very smooth. As already stressed, this cluster sample has been 
chosen for the regularity of its members. In order to asses in a significant way the effects 
of dynamical evolution on the shape of the metallicity profiles it is necessary to perform a 
statistical analysis over a large (say ^ 40) cluster sample. This task is left to a future paper, 
where a number of issues will be investigated with a statistically robust cluster sample. For 
the sake of clarity in panel (b) are also shown the profiles of the runs L40A and L40A2, to 
be compared with that of L40A3. These profiles are clearly below the measured values and 
confirm the parameters of model A3 as the ones yielding the best agreement with data. 

The iron profile of L119A3 is inconsistent with the data points; the iron abundances 
are smaller at all the radii. The profile is very similar to the ones of L40A2 and L40A. 
This is a failure of model A3 that is hard to reconcile within the framework of the adopted 
prescriptions, unless the diffusion of metals in the ICM is characterized by a minimum 
diffusion length. The most important difference of ACDM119 is that this is the coldest 
of the four test clusters, with a mass- weighted temperature of ~ IkeV. For this cluster a 
proper comparison of the simulated iron profiles with the mean metallicity profile of the 
sample is not possible. The averaged profile is that of 9 cooling flow clusters, with minimum 
temperatures ^ 4keV. Without measured profiles for cool clusters is therefore difficult to 
put observational constraints on different model parameters from the simulated profiles. 

Global values, Ape = Mpe/Mn, of the iron abundances for the simulated clusters are 
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Figure 5. Comparison against data of several observational variables at 2 = from the runs with an Arimoto-Yoshi IMF (A), 
(a): The projected emission weighted temperature profiles as a function of r/r200 for the runs A3 (Table 4); data points from 
De Grandi & Molendi (2002), Tjqq is T200 normalized to match the last data point (T200 — 8090^200 )' 0^)- panel (a), 

but for the projected metallicity profiles; data points from the CF clusters of De Grandi & Molendi (2001). Blue lines are the 
profiles corresponding to the runs of panel (a), (c): Average iron abundances at r = 0.5h~^ Mpc, open squares refer to A3 runs 
and open circles to A runs; filled squares from the nearby cluster sample of Matsumoto et al. (2000). (d): Bolometric X-ray 
luminosities as a function of the temperature. Data points from Fig. 11 of Tozzi & Norman (2001). For the sake of clarity only 
a fraction of the data set has been plotted. The continuous line is the best-fit Lx = 3.11 ■ W^"^ (T /GkeV)'^'^'^ of Markevitch 
(1998), dashed lines are the 68% confidence intervals. Filled squares are the values of Lx for the A3 runs, the filled circle refers 
to the A4 run (esN = 10^" erg). As in Markevitch, Lx for the simulated clusters has been estimated by removing a region of 
radius 50h~^ Kpc centered on the peak of the gas density. Mass- weighted temperatures have been used as estimators of the 
spectral-fit temperatures. 



compared in panel (c) for models A and A3 against the estimated values for the nearby cluster 
sample of Matsumoto et al. (2000) sample. These values have been taken from column 4 
of Table 1 of the Matsumoto et al. paper and are estimated at a radius of 0.5h~^Mpc. For 
model A3 there is a fair agreement with data, while model A yields values of Ape which are 
marginally consistent with the allowed uncertainties ( 90% c.l.). The run L119HA3 has a 
value Aj7e ~ 0.2, about a factor ~ two smaller than that expected extrapolating the sample 
average below the minimum 2keV ) cluster sample temperature. This is connected with 
the results shown in the previous panel. Without a minimum smoothing length there is a 
clear tendency in the simulations to produce a lesser amount of iron than that inferred from 
observations for cool clusters. 

The open diamond of Fig. 5(c) is the value of Ape ( ~ 0.47, see Table 5) for the run 
L119HA5. For this simulation, a minimum value /i^m — ^'^Kpc has been assumed for the 
smoothing length of metals in the ICM. Such a high value of Ape is clearly more indicated to 
account for the iron abundance of low temperature clusters. The observational evidence of an 
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Figure 6. The final profiles of the simulations L119A3 and L119A4 are compared against those of the corresponding high- 
resolution runs (thick lines), (a): radial temperature profiles, (b) : Gas density profiles, the profiles with p{0)/pc ~ 10^ are 
those of the gas metallicity. (c) : SF rates versus age. (d): Iron abundance profiles. 



iron abundance decreasing with Tx is however statistically weak (Mushotzky & Lowenstein 
1997; Finoguenov, Arnaud & David 2001). Furthermore, an increase of Ape for cool clusters 
can be an artifact due to the presence of a dominant galaxy in the cluster central region. For 
example, Fukazawa et al. (2000) have removed the contribution of the central region for those 
clusters with a cD galaxy and for their cluster sample found that there is not a significant 
correlation Ap^. — kTx- This is also in agreement with what was found by Finoguenov et al. 
(2001). A statistical comparison between the Ap^. distribution generated by the simulations 
and the one from real data has been performed for the models discussed above. The statistical 
tests applied to the Ape data and the Ape distribution of the simulated cluster sample are 
the Student t-test for the means, the F-test for the variances and the Kolgomorov-Smirnov 
(KS) test for the distributions. The corresponding probabilities give the significance level 
that, according to the analyzed quantities, the two sets are originated from the same process. 
The results are presented in Table 6. Model A is clearly ruled out, model A3 is marginally 
consistent. This is because the value Ap^, of L119HA3 is significative to lower the mean of 
the sample. A much better agreement is obtained if the global iron abundance is given by 
the run L119HA5 for the cluster ACDM119. 

Finally, the final values of the bolometric X-ray luminosity are shown in panel (d) as a 
function of the cluster temperature. Mass-weighted temperature have been used as unbiased 
estimators of the spectral temperatures ( Mathiesen & Evrard 2001). The luminosity Lx of 
the simulated clusters is calculated as described in sect. 2. Data points are those of Fig. 11 
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data-model 


Pt 


Pf 


Pks 


Ap^-A 


.003 


.032 


.001 




.082 


.694 


.120 




.689 


.993 


.756 


/9as(500)-A3 


.02 


.565 


.018 


fgas{200)-A3 


.001 


.221 


.001 


/sa4500)-A3(19)'' 


.143 


0.849 


0.163 



Table 6. Statistical tests applied to distributions obtained from real data and simulated clusters, pt is the 
Student t-test applied to the means, pp is the F— test for variances and Pks is the KS statistic to discriminate 
two distributions. Ap^ is the measured iron abundances distribution of Fig. 5, panel (c). fgas refers to the 
observational sample of the two top panels of Fig. 8. The parameters of the model are indicated in Table 4. 
For the simulated cluster sample the numerical resolution of the runs is LOO, L39, L40 and L119H. This 
sample has L119HA5 in place of L119HA3. ^: This sample has only three clusters, but their initial conditions 
have been set with Sl[,h^ = 0.019. 

of Tozzi & Norman (2001). For the sake of clarity, not all the points of the Figure are plotted 
in the panel. For a consistent comparison with data, a central region of size 50h~^Kpc has 
been excised (Markevitch 1998) in order to remove the contribution to Lx of the cooling 
flow central region. For model A3, the Lx of the simulations are in excellent agreement with 
data over the entire range of temperatures. An additional run has been performed for the 
cluster with the lowest temperature (ACDM119). The parameters of this run (A4) are the 
same of model A3, expect for a SN feedback energy of 10^° erg being used for both SNII and 
la. This value is 10% of that of model A3 and has been considered in order to investigate 
the effects on final X-ray properties of the amount of heating returned to the ICM by SNe. 
As can be seen, for model A4 the final Lx of the run L119HA4 is very similar to that of 
L119HA3. The result demonstrates that final X-ray luminosities of the simulations are not 
sensitive to the amount of SN feedback energy that has heated the ICM. In fact, a simulation 
run with a zero SN energy returned to the ICM {esn = 0) yields very similar results. These 
results are particularly relevant in connection with the recent proposal (Bryan 2000) that 
the X-ray properties of the ICM are driven by the efficiency of galaxy formation, rather 
than by heating due to non-gravitational processes. Final profiles for models A3 and A4 are 
investigated in more detail in Fig. 6 and 7. 

4.4 SN heating of the ICM and cooled gas fractions 

In Fig. 6, the final profiles of the runs L119A3 and L119A4 are compared against those 
of L119HA3 and L119HA4. The main differences are between the profiles of the standard 
resolution run and the high-resolution run. The profiles of the runs for the models A3 and 
A4 are very similar, with the only important difference between the two models being the 
SFR at early redshifts. The temperature profiles are almost identical, so energy feedback 
from SNe is not relevant to determine final gas properties. This conclusion is also supported 
from the entropy profiles of Fig. 7. 

In Fig. 7, final distributions of the two runs L119HA3 and L119HA4 are compared. In 
panel (a) are plotted the Si abundances synthesized in the explosions of SNII and SNIa. 
There are not large differences between the corresponding profiles. The Si abundances of 
SNIa are almost identical. For SNII, the Si abundances profile of L119HA4 is higher than 
that of L119HA3 in the cluster central region ( ^ O.lMpc). These differences follow because 
the run L119HA4 has at early redshifts a higher SFR than that of L119HA3 (see also Fig. 
6c). This is a consequence of the reduced thermal pressure and higher gas densities in the 
former simulation. This effect is not relevant for the SNIa explosions because the involved 
timescales are much higher than those of the early generations of SNII. 
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The cumulative injected energy per baryon within the radius r is defined as the ratio 
Esn{< r)fi'mp/Mg{< r), where Esn{< is the total SN energy injected within that radius, 
Mg{< r) is the corresponding gas mass. This ratio is shown in panel (d) for the two runs. 
The ratio has a radially decreasing behavior, for the run L119HA3 drops from ~ IkeV/part 
at r ~ O.lMpc to ~ 0.5keV/part at r 2^ IMpc ~ r2oo- There is a small drop in the 
cluster inner region, presumably due to a biased estimate because of the small number of 
simulation particles for r ^ AOKpc. For the run L119HA4 the energy per particle ratio has 
a shape similar to that of the corresponding ratio for the run L119HA3, at r ~ IMpc takes 
the value of ~ O.OAkeV/part, which is about ~ 10% of the value of the parent run with 
EsN = lO^^erg. The binned distribution (panel c) is much more noisy than the cumulative 
one, but it qualitatively confirms the expected behavior. 

In order to assess the amount of heating from the SNe, this ratio is often inferred from 
the measured abundance of metals (Finoguenov, Arnaud & David 2001). Si abundances are 
particularly useful because of the similar yields for different SN types. The ratio estimated 
from the amount of Si abundance Asi is esN-AsiUmp/ysi with ysi ~ 0.12 being the average 
SN yield of Si. The estimated ratio (thin lines) is in good agreement at all the radii with the 
values obtained from the simulation. For Esn = lO^^ergf the average energy per particle of 
the cluster ACDM119 is then ~ {^.bkeV/part at the virial radius (r2oo ~ IMpc). This cluster 
has a virial temperature ~ l.bkeV, from a sample of 18 relaxed clusters with temperature 
below ~ AkeV. Finoguenov, Arnaud & David (2001, see Fig. 9) find an average SN injected 
energy of ~ 0.5keV/part for a cluster with a temperature ~ IkeV. There is therefore a good 
agreement between the measured amount of thermal energy per particle associated with SN 
feedback and that predicted by the simulations with the model parameters of Fig. 7. 

The gas entropy is defined as S{r) = kBT{r)/ne{ry^^, where rig is the electron number 
density. For the two runs L119HA3 and L119HA4 the entropy profiles in units of keVcm'^ are 
plotted in panel (b) of Fig. 7. An important result that follows from the radial dependence 
of the two profiles is that they are almost identical. A simulation with a zero SN energy 
yields a very similar profile. This means that final ICM properties are largely unaffected by 
the amount of energy injected by SNe into the ICM. This follows because most of the energy 
is injected in the cluster central region where the gas density is higher and cooling is very 
efficient to radiate away the energy of the reheated gas. Energy feedback from SNe can modify 
the ICM state for small clusters or groups with Tx ~ IkeV. Another important result from 
the entropy profiles of Fig. 7 is the entropy level of the gas at r ~ O.lRv — 0.1r2oo — O.lMpc. 
At this radial distance the gas entropy of the simulated cluster is ~ 200keVc'm'^, a value 
consistent with a set of observations (Ponman, Cannon & Navarro 1999; Lloyd-Davies, 
Ponman & Cannon 2000; Wu & Xue 2002a) for a system with Tx ~ 1.5keV. Together with 
the good agreement between the Lx — Tx relation obtained from simulations with that from 
observations these findings provide strong support for the radiative cooling model proposed 
by Bryan (2000). 

According to this scenario (Bryan 2000; Voit & Bryan 2001; Wu & Xue 2002a; Muanwong 
et al. 2002 ; Dave et al. 2002) the efficiency of galaxy formation is higher in groups and 
cool clusters than in hot clusters. The low-entropy cooled gas is removed because of galaxy 
formation, and is replaced by the infiow of the surrounding high-entropy gas. The higher 
efficiency of galaxy formation for cool systems explains the central entropy excess over the 
self-similar predictions (Ponman, Cannon & Navarro 1999). One of the most important 
observational consequences of the cooling model is that the fraction of hot gas fg = Mg/M^ 
increases going from cool clusters to hot clusters. Conversely, the fraction fstm- of cooled gas 
turned into stars should decrease with the mass of the system. 

Observational evidence for a dependence of fg with Tx is controversial (David et al. 1990; 
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Figure 7. For the cluster ACDM119 final distributions of the run L119HA3 are compared against those of L119HA4 (esN = 
10^'^ erg). Panel (a) shows the Si abundance of the gas, originated from both SNII and SNIa metal enrichment of the ICM. 
(b): Entropy profiles S = T{r)/ne{r)^/^ , where rxe = pg{r)/ fie^ip and /ie — 1.14. (c): Distribution of the gas particle energy 
injected by SNe. (d): As in panel (c), but for the cumulative distribution -Esiv(< ^)/^p(< where Np = Mg{< r)/fj,mp. 
Thin lines correspond to Egpf{< r)/Np{< r) as obtained from Si abundance. 

David, Jones & Forman; Mohr, Mathiesen & Evrard 1999; Arnaud & Evrard 1999; Roussel, 
Sadat & Blanchard 2000; Balogh et al. 2001). The main difficulty is the extrapolation to 
virial radii of the X-ray data for cool clusters, which can bias the estimates of the cooled 
gas fraction. According to Roussell et al. (2000) there is not strong observational support 
for an increase of fgas with Tx- This is in disagreement with the findings of David, Jones & 
Forman (1995), who show that the gas fraction fgas increases with the X-ray temperature. 
According to Mohr, Mathiesen & Evrard (1999) there is a weak dependence of fgas with Tx 
in their X-ray ROSAT sample of 45 galaxy clusters. A constant fgas is inconsistent with the 
data at 95% confidence level. 

Arnoud & Evrard (1999) have analyzed the temperature dependency of fgas for a sam- 
ple of 24 clusters with accurate measured temperatures and low cooling flows. They have 
determined fgas for two radii enclosing gas overdensities S = 500 and 6 = 200. Total masses 
have been estimated according to two different methods. The /3-model (BM) assumes an 
isothermal gas density proffie in hydrostatic equilibrium with the shape determined accord- 
ing to the X-ray surface brightness. The virial theorem (VT) estimates the cluster total mass 
according to virial equilibrium at a fixed density contrast. The relation is calibrated from a 
set of numerical experiments (Evrard, Metzler & Navarro 1996) and is not sensitive to the 
assumed cosmological model. According to Arnaud & Evrard (1999) it is difficult to draw 
conclusions on the dependence of fgas on the cluster temperature. The results depend on 
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the chosen model. For the BM model, the difference in the value of fgas between subsamples 
of cool clusters and hot (> AkeV) clusters is not statistically significant. For the VT model 
this difference is at the Scr level at a radius corresponding to a gas overdensity 5 = 500. 
Arnaud & Evrard (1999) indicate a conservative upper limit of ~ 30% in the la fractional 
error of fgas- 

For the cluster sample of Arnaud & Evrard (1999), the gas fraction fgas is shown as 
a function of the total cluster mass in the top panels of Fig. 8. The sample distributions 
are derived from the VT model at a density contrast of 6 = 500 (top left panel) and 
6 = 200 (top right panel). The data points are those in the top panels of Fig. 3 of Arnaud 
& Evrard (1999), with an assumed ^ 30% uncertainty, (top right panel). For the model A3 
the corresponding values of fgas at 5 = 500 and S = 200 from the numerical cluster sample 
are also shown as open squares. There is a clear tendency for the two simulated cluster 
distributions to follow the fgas distribution of the data, but with a reduced amplitude. The 
result of a statistical comparison are reported in Table 6. The fgas distribution of model A3 at 
6 = 500 is inconsistent with data with a high significance level (95% c.l.), and the situation at 
6 = 200 is even worse. For this overdensity, an extrapolation of X-ray data up to the required 
radius has been performed (Arnaud & Evrard 1999) in order to evaluate the corresponding 
data points. Therefore, possible biases can undermine the estimate of fgas at 5 = 200. For 
fgas[S = 500] a possible source of disagreement between the numerical distribution and 
that of data points lies in the assumed value of the cosmological baryonic density Qb in 
the simulations. Here it has been assumed Q^h"^ = 0.015, but recent measurements (Buries 
& Tytler 1998) favor Q^h"^ = 0.019. In correspondence to this value of Qf, in the initial 
conditions of the simulations, the values of fgas[S = 500] are displayed in Fig. 8(a) for a 
subsample of three simulated cluster (ACDMOO, ACDM119 and a new cluster with virial 
mass ~ 2.5 ■ lO^^/i~^M0). The distribution is now in a better agreement with the data, and 
the confidence levels for rejecting the null hypothesis are now below 95%. Therefore it is fair 
to say that the ICM gas fractions obtained here are consistent with available observational 
estimates. 

The predictions of the radiative cooling model for the ICM evolution have been recently 
investigated in a number of papers, either through analytical methods (Bryan 2000; Voit 
& Bryan 2001; Wu & Xue 2002a; Voit et al. 2002; Wu & Xue 2002b) or with numerical 
simulations (Muanwong et al. 2001; Muanwong et al. 2002 ; Dave et al. 2002). Results 
from the employed methods clearly show a decreasing fstar with the cluster temperature. 
The dependence of fstar[S = 200] on Tx is shown in Fig. 8(c) for the four test clusters. 
Open squares are for model A3 with standard resolution and filled squares refers to the 
high-resolution runs. For comparative purposes the fstar predicted by the analytical model 
of Bryan (2000) is also shown. Because of the different values of Qb and Qm in the paper of 
Bryan, fstar has been approximately rescaled by the factor fbihere) / fb{Bryan) = 0.1/0.16 = 
0.625. There is a fair agreement of the model with the simulations at high temperatures, 
but for Tx ~ IkeV fstar of the simulations is below the model predictions by a factor ~ 
two. The disagreement at low temperatures is ameliorated for the runs with Qbh"^ = 0.019, 
in such a case fstar{Bryan) must be rescaled by the factor ~ 0.81. The predicted value of 
fstar is now ~ 0.075 at Tx — IkeV, which is a factor ~ 50% higher than the value ~ 0.05 of 
fstar at the same temperature for the runs with Qbh"^ = 0.019. 

The comparison of fstar with data is a controversial issue (Balogh et al. 2001), also 
because of the lack of firm bounds on the allowed uncertainties (Roussel, Sadat & Blanchard 
2000); however, it is worth stressing that Roussel et al. (2000) found a weak dependence of 
the gas fraction on cluster temperature when the estimates of virial masses are calibrated 
from numerical simulations rather than inferred from the /3-model. The dependence is not 
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Figure 8. Top panels: gas fraction fg = A/g(< r)/MT{< r) versus Mt within the radius within which p/pc = S, (a) is for 
S = 500 and (b) is for S = 200. Data points from Fig. 3 of Arnaud & Evrard (1999, VT model). The open circles in panel 
(a) are the values of fg for runs with fli, = 0.019h~^ . (c) : The ratio of the mass of cold gas turned into stars to the total 
mass is shown as a function of the cluster temperature. The dashed line is the Bryan (2000) best-fit fstar to data, rescaled by 
the factor fbaryon{here) / fiiaryon{Bryan) = 0.1/0.16. (d): Prom the left, gas mass above a given temperature T is plotted as 
a function of the temperature itself. Dashed lines axe for the very-high resolution runs VH, while L119HA3 is for ACDM119. 
From the right; the fraction of total X-ray luminosity that is emitted from the gas below a given temperature is shown for the 
runs corresponding to the gas masses on the left. 



statistically significant because of the large size of the sample error bars. A comparison of 
fstar predicted by an analytical model has been performed by Wu & Xue (2002b) against 
the Roussell et al. (2000) data. The theoretical predictions are of about twice as high as the 
sample values. From Fig. 8(c) the values of fstar[S = 200] for the three runs with Qbh"^ = 0.019 
can be compared with the analytical predictions of the Wu & Xue (2002b) model (cf. Fig. 2 
of their paper). The values of Qf, and h are slightly different, but there is a broad agreement. 
At Tx ~ IkeV fstarihere) ^ 0.05 is lower than the theoretically predicted value ~ 0.07. As 
noted by Wu & Xue theoretical models of cooling yield a dependence of fstar on Tx steeper 
than that found in hydrodynamical simulations and observations. At low temperatures fstar 
of the simulations is in better agreement with the Rousell at al. (2000) data, but the observed 
uncertainties do not allow firm conclusions. The simulations results for fstar are also in 
qualitative agreement with those of the radiative model of Muanwong et al. (2002). In units 
of l^fe/l^m, the values of fstarifg) range from ~ 0.38(0.46) from M200 = 0.8 ■ IO^^/i-^Mq to 
~ 0.34(0.56) for a cluster with a virial mass of M200 ~ 2.4- 10^^h~^MQ. These values are not 
very different from those which can be inferred from the distributions plotted in Fig. 3(b) 
on Muanwong et al. (2002) for their radiative model. The cosmological parameters of the 
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simulations of Muanwong et al. (2002) are identical to those assumed here {Vt^ is a ~ 20% 
higher), and the numerical resolution is broadly similar. 

Finally, the effects of the numerical resolution on the predicted final amount of cooled 
gas are shown in panel (d) of Fig. 8. As a function of the gas temperature the fractional 
contribution to the total X-ray luminosity is plotted from the right for runs of three clusters 
(ACDMOO, ACDM40 and ACDM119) with different resolutions. For the first two clusters the 
ratio of Lx{< T)/Lt for very high-resolution runs are compared against that of standard 
resolution. The ratio Lx{< T)/Lt of the run L119A3 is instead compared against that of 
L119HA3. From the left, the gas masses above a given gas temperature are plotted for the 
corresponding runs. The distribution of the gas mass versus the gas temperature shows that 
resolution effects are most important in the low-temperature region of the gas temperature 
distribution, whereas X-ray luminosities are largely determined by the high-temperature 
part of the distribution. This demonstrates that the fraction of cold gas depends on the 
numerical resolution, as can be seen from Fig. 8(c), but is not a source of concern as far as 
X-ray properties are interested. 

For the cluster ACDM119, the quantity, Mgas{> T) is weakly dependent on the numerical 
resolution of the simulations. From Fig. 8(c), one can see that fstar at low temperatures also 
seems to converge as the resolution is increased. This issue is strictly related to the global 
value of the baryonic fraction of cooled gas Vtcoid/^b = fbid^obal). Observational upper 
limits indicate f},{global) < 10% (Balogh et al. 2001). According to Balogh et al. (2001), the 
global fraction fi,{global) is expected to increase as the resolution limit of the simulation is 
increased, unless a feedback model is incorporated in the simulation. The results of the high- 
resolution runs suggest that at low temperatures convergence is being achieved for fstar and 
that the feedback model implemented here effectively regulates the amount of cooled gas. 
This is strongly indicated by the large value of fstar (~ 0.05) when the SN feedback energy 
is reduced (open diamond of Fig. 8(c), which corresponds to the run L119HA4). For clusters 
with temperatures above ^ 2keV, the values of fstar appear to converge for simulations with 
a number of gas particles Ng ^ 200,000. This is shown in Fig. 8(c), where for the cluster 
ACDMOO the value of fstar ioT the very-high resolution run L00VHA3 (filled triangle) is 
~ 0.03 and is close to that of fstar for L00HA3. Note that the mass of the gas particle for 
the run L00VHA3 is only a ~ 10% higher than that of L119HA3 (~ 7 ■ IO^Mq). These 
issues can be clarified in deeper details with a set of cosmological simulations with different 
resolutions that incorporate the SN feedback scheme implemented here. Such analysis is 
beyond the scope of this paper and will be considered in a future work. 

5 SUMMARY AND CONCLUSIONS 

In this paper results from a large set of hydrodynamical SPH simulations of galaxy clusters 
have been used to investigate the dependence of iron abundances and heating of the ICM 
on a number of model parameters. The simulations have been performed with different 
numerical resolutions and a numerical cluster sample covering nearly a decade in cluster 
mass. The modeling of the gas physical processes in the simulations incorporates radiative 
cooling, star formation and energy feedback. The gas is metal enriched from SN ejecta of 
type II and la, and the cooling rate depends on the local gas metallicity. This allows us to 
follow self-consistently the ICM evolution also for cool clusters, whose departures of cluster 
scaling relations from self-similarity are most important. 

The metal enrichment of the ICM is governed by a number of model parameters. Theo- 
retical uncertainties on the shape of the IMF and nucleosynthesis stellar yields lead to final 
iron abundances which can differ by a factor ~ two. These two parameters have then been 
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kept fixed for a large fraction of the performed simulations. This is done in order to restrict 
the full range of the parameter space. They have been chosen so that the corresponding 
final amounts of iron in the ICM are the ones that most closely agree with that indicated 
by observations. The ejected metals are distributed among gas neighbors according to the 
SPH formalism. Final iron abundances have been found to be sensitive to the shape of the 
smoothing kernel of the ejected metals Wz and to the constraints on the corresponding 
smoothing lengths h'^ . For a certain choice of the above parameters, it has been found a 
fair agreement between the observed radial metallicity profiles with those obtained from the 
simulated clusters. Global iron abundances are also consistent with measured values. There 
are however a number of issues which are still open. 

i) For cool clusters the ejection parameters are not well determined because of the lack 
of measured iron profiles, also constraints from global iron abundances suffer from selection 
biases in the cluster samples. The simulated profiles have a radial decay steeper than ob- 
served. This difference could be mainly due to the lack in the simulated cluster sample of 
mergers which have efficiently remixed the metal content of the ICM. The small size of the 
numerical cluster sample does not allow to reach firm conclusions. In order to investigate 
the correlation between the shape of the final metallicity profiles with the cluster dynamical 
evolution it is then necessary to analyze simulation results from a statistically robust (say 

^ 40) cluster sample. 

ii) A discrepancy with observations is given by the radial behavior of the projected 
emission-weighted temperature profiles, which have a steep rise at r ^ 0.1r2oo- This feature 
is not shared by observations, for the innermost bin the simulated temperatures are higher 
by a factor ~ two than the measured values. This difference cannot be easily accommodated 
by the gas physical modeling of the simulations. Because of radiative cooling, an increase 
of the gas temperature toward the cluster centre is expected as a consequence of entropy 
conservation (Wu & Xue 2002a). A certain amount of cooled gas is present in the core of 
the clusters and is responsible of the strong decline of the temperature profiles at distances 
very close to the cluster centres (r ^ 0.02r2oo)- This feature could help to reconcile the 
differences with the shapes of the observed profiles. As found by Mathiesen and Evrard 
(2001), the measured spectral fit temperatures can be biased toward lower values because 
of line cooling. An X-ray spectral fit analysis is then necessary to resolve the discrepancies 
between the simulated and measured temperatures profiles in the cluster core regions. 

The ICM X-ray properties of cool clusters are largely unaffected by the amount of feed- 
back energy injected by SNe, with the luminosity-temperature relation being in excellent 
agreement with data. The final gas entropy distribution has been found almost independent 
on the heating efficiency, for Tx ~ IkeV the core entropy excess is in agreement with es- 
timated values. These findings support the radiative cooling model of Bryan (2000), where 
the ICM X-ray scaling relations are driven by the efficiency of galaxy formation. The model 
predicts that at the cluster virial radius the fraction of hot gas fg is positively correlated with 
the cluster temperature Tx- A number of observational uncertainties prevents to put tight 
constraints, but the fg — Tx distribution of the simulated cluster sample is not inconsistent 
with available estimates from X-ray data. 

To explain the steepness of the luminosity-temperature relation an alternative to the 
cooling model is the preheating scenario, where the gas has been heated by an energy 
injection that occurred at early epochs. Simulation results (Bialek, Evrard & Mohr 2001; 
Muanwong et al. 2002 ; Borgani et al. 2002) show that a number of observed X-ray properties 
can be reproduced in the preheating model. Future X-ray observations of cool clusters and 
groups with Chandra will help to discriminate between the two models. It is worth noting 
that simulations of galaxy clusters in the preheating scenario must also include radiative 
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cooling and SN feedbacks. Therefore, the agreement with observations of X-ray properties 
reproduced here by the simulated clusters in the cooling model suggests that additional heat 
sources, if required, are most likely to affect the ICM properties of systems with temperatures 
below ^ IkeV. 
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